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Abstract 

If dark matter particles pair annihilate into stable Standard Model 
particles, a population of energetic, non-thermal electrons and positrons 
is necessarily produced. The secondary radiation resulting from the 
energy losses of this non-standard population (including synchrotron 
and inverse Compton up-scattering of background radiation photons) 
in turn produces a peculiar multi-wavelength spectrum, extending 
from radio to gamma-ray frequencies. We give here an overview of 
multi-wavelength searches for dark matter, including the expected in- 
jection spectrum and production rate of electrons and positrons from 
dark matter annihilation, the computation of the effects of propaga- 
tion and energy losses and the actual multi-wavelength emission. We 
then outline the application of this general framework to the case of 
galaxy clusters (specifically Coma, IE 0657-56 - the so-called Bullet 
cluster - and Ophiuchus) and of dwarf spheroidal galaxies (including 
Draco, Fornax, Ursa Minor and Carina). We also review the applica- 
tion of multi-wavelength search strategies to our own Milky Way, and 
more specifically to the Galactic center environment, to dark galactic 
mini-halos and to the so-called WMAP haze and other radio data. We 
argue that multi-wavelength observations will complement gamma-ray 
observations as probes of particle dark matter, since the expected lu- 
minosities at different frequencies are generically comparable. The 
indirect search for dark matter with astronomical observations at var- 
ious frequencies is therefore a crucial and indispensable element in the 
quest for the fundamental nature of dark matter. 

^ This contribution appeared as chapter 27, pp. 547-564, of 'Particle Dark 
Matter: Observations, Models and Searches' edited by Gianfranco Bertone, 
Copyright 2010 Cambridge University Press. Hardback ISBN 9780521763684, 
http://cambridge.org/us/catalogue/catalogue. asp ?isbn=9780521 763684 
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1 Introduction 



The class of Weakly interacting massive particles (WIMPs) is the lead- 
ing category of particle dark matter candidates: On the one hand, 
the mechanism of thermal freeze-out of a stable WIMP x leads to 
a non-relativistic relic population whose relative matter density can 
be approximated as ^i-^/i^ ~ 3 x 10^^'^ cm^s~"'^/((JAf), where h is the 
Hubble constant in units of 100 km s~^ Mpc~^, fi^ is the ratio of the 
X density over the critical density, and {(tav) is the thermally aver- 
aged X pair-annihilation cross section. If the new physics connected 
to the WIMP is at the electro- weak scale, msW) one can estimate 
(aAv) w ckew/^^ew ~ 10^^^ cm^s^^, giving a thermal leftover x 
population that lies in the same ballpark as the dark matter abun- 
dance inferred from cosmic microwave background anisotropics, large 
scale structure and other astronomical observations |22j ; On the other 
hand, numerous motivated particle physics extensions to the Standard 
Model encompass a stable WIMP, including supersymmetry [2] and 
universal extra-dimensions (UED) [16j, in virtue of unbroken discrete 
symmetries (i?-parity in the case of supersymmetry, and Kaluza-Klein 
parity in the case of UED). 

Since WIMPs were once kept in thermal equilibrium by pair anni- 
hilations into Standard Model particles and inverse WIMP pair pro- 
duction processes, even in today's cold Universe, occasionally, WIMPs 
can pair annihilate, giving rise to energetic, stable "ordinary" Stan- 
dard Model particles. The rate of WIMP pair annihilation depends 
on the number density of dark matter in a given location and on the 
pair-annihilation rate {aAv)- The spectrum and particle content of 
the final state products of a pair-annihilation event are, instead, de- 
termined once the microscopic details of the specific WIMP model 
are specified. Given a WIMP x, the production rate of a final state 
Standard Model particle species / = 7, e^,i?, v, ■ ■ emitted in a pair- 
annihilation event in a location x with dark matter density /Odm(^)) 
reads 

<3(j..*) = 4^(...»>fi. (1) 

In the formula above, dA'^^ is the number of / particles with modulus 
of the momenta between p and p + dp resulting from the WIMP pair- 
annihilation event. Numerous studies have investigated the indirect 
detection of WIMPs through signatures in the spectrum of cosmic-ray 
positrons and antiprotons, or through the emission of energetic gamma 
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rays. With the recent successful launch of the space-based antimatter 
detector PAMELA [T7] and of the Fermi gamma-ray space telescopes 
|19| . dark matter indirect detection faces new opportunities, with a 
forthcoming dramatic increase in the quality and quantity of available 
data. 

The pair-annihilation of WIMPs yields gamma rays from the two 
photon decay of neutral pions, and, in the high energy end of the spec- 
trum, from internal bremsstrahlung from charged particle final states. 
Sub-dominant one-loop annihilation modes give rise to the additional 
emission of monochromatic photons: the direct pair-annihilation into 
two photons, for instance, yields two gamma rays with = m^. 
Photons from prompt emission are generally in the gamma-ray band, 
where absorption is negligible, hence the induced fluxes or intensities 
can be straightforwardly derived summing contributions along the line 
of sight. 

Concerning lower frequencies, a conspicuous non-thermal popula- 
tion of energetic electrons and positrons results from the decays of 
charged pions produced by the hadronization of strongly interacting 
particles in the final state, as well as from the decays of gauge bosons, 
Higgs bosons and charged leptons. This non-thermal population 
looses energy and produces secondary radiation through several pro- 
cesses: synchrotron in the presence of magnetic fields, inverse Comp- 
ton scattering off starlight and cosmic microwave background photons 
and bremsstrahlung in the presence of ionized gas. This radiation can 
actually cover the whole electromagnetic spectrum between the radio 
band to the gamma-ray band. 

The computation of the multiwavelength emissions from WIMP- 
induced energetic involves several steps: (1) the assessment of the 
spectrum and production rate of electrons and positrons from dark 
matter pair annihilations (that follows from Eq. ([T| with / = e^); 
(2) the computation of the effects of propagation and energy losses, 
possibly leading to a steady state configuration; (3) the compu- 
tation of the actual emissions from the mentioned equilibrium con- 
figuration; (4) the evaluation of eventual absorption of the emitted 
radiation along the line of sight to derive fluxes and intensities for a 
local observer. We describe items (1) and (2) in the next section [LT| 
while we devote sec. 11.21 to the calculation of the secondary emission 
from dark matter annihilation induced e^. Sec. [2] then applies the 
general outlined framework to the case of galaxy clusters, speciflcally 
Coma (sec. [2T|, IE 0657-56 (the so-called Bullet cluster, [22]) and 
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Ophiuchus (|2.3p . The following sec. [3] studies the other extreme of 
galactic scales: dwarf galaxies. In particular, we focus on the case of 
the dwarf spheroidal (dSph) galaxy Draco in l3.1l and on constraints on 
dark matter particle models from X-ray observations of dSph galaxies 
Fornax, Ursa Minor and Carina in 13.21 Sec. U] specializes to the case of 
our own Milky Way, and more specifically on the Galactic center (j4.ip . 
on dark galactic mini-halos (j4.2p and on the dark matter annihilation 
origin of the so-called WMAP haze ()4.3p . Finally, sec. [5] addresses 
radio constraints from cosmological dark matter annihilation at all 
redshifts, and sec. [6] concludes and provides an overview. 

1.1 The electron-positron equilibrium number 
densities 

Except for WIMPs which annihilate mainly into light leptons, the bulk 
of electrons and positrons from WIMP annihilations are secondaries 
from intermediate steps in which charged pions are produced and de- 
cay along the chain: vr^ — )■ ^^v^{i>^j), with — )• + ^^(t'^) + fe(^'e)- 
These are twin processes with respect to the production of neutral 
pions and their vr'^ — )• 77 decay, hence at the level of the production 
rates, the yield into electron and positrons in Eq. [1] is at the same 
level of the gamma-ray yield and with analogous spectral shape. In 
principle one should take into account other stable yields, such as 
protons, which can be confined in the astrophysical environment in 
which WIMP annihilation is considered, and subsequently give rise to 
further secondaries (such as pions in proton collisions with the ambi- 
ent interstellar material); their multiplicity in WIMP annihilation is 
however much smaller than for electrons and positrons and, therefore, 
they can be safely neglected. 

After emission, electrons and positrons go through a random walk 
in the turbulent and regular components of ambient magnetic field, 
loosing energy by radiative emissions, as well as loosing and/or gain- 
ing energy in the convective or advective effects associated to plasma 
outflows or inflows. The sum of these effects are usually described in 
terms of a transport equation of the form: 





dt 



p 



(2) 
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where ne{r,p,t) is electron/positron number density per unit parti- 
cle momentum p, Dxx is the spatial diffusion coefficient, diffusive re- 
acceleration is described as diffusion in momentum space and is de- 
termined by the coefficient -Dpp, p = dp/dt is the momentum gain or 
loss rate, and V is the convection or advection velocity. Except for 
when considering the effect of individual dark matter substructures 
orbiting within a given astrophysical object the WIMP source func- 
tion can be regarded as time independent, and the emission rate is 
large enough to reach equilibrium, so that Eq. [2] can be solved in the 
steady state approximation, i.e. setting the left-hand-side of Eq. [2] to 
zero. Simplifying ansazte on the morphology of the source (namely 
spherical symmetry when considering galaxy clusters, dwarf satellites 
and the galactic center region, axial symmetry when considering the 
whole Milky Way) and on the spatial dependence of the parameters 
in the model are usually adopted; another useful approximation is to 
consider the limit of free-escape at the boundary of the diffusion re- 
gion, i.e. impose rig = outside the region where it is estimated that 
the propagation equation can be applied. We discuss briefly below 
how relevant of the different terms in the transport equation are for 
astrophysical sources we will consider. 

1.1.1 Energy loss term 

The radiative losses affecting the e"^ — e~ propagation are inverse 
Compton scattering on CMB and starlight, synchrotron emission, 
bremsstrahlung, ionization and Coulomb scattering. The term p in 
Eq. [2] can be written as: 

P = PIC (P) + Psyn {p) + Pbrem {p) + Pion {p) + PCoul {p) 
- Picip) +Psyn{p) 

where we have indicated that in all cases of our interest, the relativis- 
tic regime, either inverse Compton or synchrotron dominate, and the 
coefficients b^^ ~ 0.76 and b^y^ ~ 0.025 are found, respectively, for 
a background photon energy density Uph of 1 eV cm^^ (the CMB en- 
ergy density is about 0.25eVcm~^, while the starlight energy density 
in the Milky Way is about 0.6eVcm~^ locally and 8eVcm~^ in the 
galactic center region [23j ) and a radiating magnetic field i? of a 1 //G 
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strength (say, within one order of magnitude of what is expected in 
galaxy clusters, dwarf galaxies or the local portion of the Milky Way, 
while in the Galactic center region magnetic fields should be much 
stronger, at the 1 G level or even higher). 

1.1.2 Spatial diffusion and re-acceleration 

There is no definite prediction for the diffusion coefficient Dxx{p)- 
Its spectral shape a and normalization depend on unknown variables 
describing turbulence, namely the amplitude of the random magnetic 
field and the scale and the spectrum of the turbulence. E.g., one 
finds that in case of Bohm diffusion a = 1, while for Kolmogorov and 
Kraichnan diffusion, the scaling in energy is, respectively, a = 1/3 
and a = 1/2. On top of this, there is the uncertainty on whether and 
how to include the re-acceleration term Dpp{p), which is in general 
predicted to be related to D^xip)- When considering the case for 
diffusion of cosmic rays within the Milky Way one can rely on the 
comparison between the model and a rather vast and variegated set 
of observables, such as the ratio of local measurements of abundances 
of secondary versus primary cosmic-ray species, and global maps of 
photon emissivities in different energy bands. A nearly self-consistent 
picture (see some more discussion below) can be obtained for, say, a 
diffusion coefficient at 1 GeV of the order of ~ lO^^cm^s"^, a ~ 0.6 
and no re-acceleration, or q = 1/3 and including re-acceleration (see. 



Assuming, for simplicity, that re-acceleration can be neglected, 
to have a feeling on whether diffusion is relevant, we can guess, on 
dimensional grounds, that the diffusion length scale should be of the 
order of the square root of the diffusion coefficient times the time scale 
relevant in the transport equation, i.e. Xxx ~ V^xxTIoss, with tioss 
the time scale for the energy loss associated to radiative processes, 
namely tioss ~ p/p- This guess gives a result that is not too far from 
the diffusion length scale one obtains in one of the limits in which 
the transport equation can be solved analytically, namely for spatially 
constant diffusion coefficient and energy loss term, and neglecting re- 
acceleration and convection [1]: Xxx{p) = '^\/v{p) — v{p'), where we 
introduced the variable v, in place of the momentum p of the radiating 
electron/positron, through the double change of variables: 



e.g. m)- 




and 
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Figure 1: The figure shows the distance (A^;)^/^ which, on average, a electron 
covers while losing energy from its energy at emission E' and the energy when 
it interacts for a few values of E: 30 GeV, 10 GeV, 5 GeV and 1 GeV, and 
for a few values of the magnetic field (in /iG); we are focusing on a WIMP 
of mass 100 GeV, hence cutting E' < 100 GeV. 

and p' is the positron/electron momentum at emission. As a sample 
case, we plot in Fig. [T]the square root of Av = v{p) — v{p'), for a few 
values of radiating and emission energies, for a few sample values of the 
ambient magnetic field in units of fiG and with diffusion coefficient 
within an extrapolation which should apply to galaxy clusters (see e.g. 
[31 [5]). Length scales of few (tens of) kpc should be compared to the 
typically scales over which the dark matter density profiles varies or 
to the length scale corresponding, for a given object distance, to the 
angular resolution of the detector searching for the induced signal. 
In galaxy clusters the dark matter halo scale parameter is generally 
much larger, and diffusion turns out to be relevant only for dark matter 
profiles which are singular towards their center, on angular portions 
on the sky which are however small compared detector resolutions; it 
follows that spatial diffusion is a marginal effect in galaxy clusters, 
and can be safely neglected for most practical purposes in this case 
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In the limit of negligible spatial diffusion (as well as negligible 
re-acceleration and convection), the equilibrium number density Ue 
reduces to: 

1 



ne{x,p) = J dp' Qe{x,p') . (5) 

On the other hand, in case of Milky Way dwarf satellites (whose size 
is few kpc), and for estimates of the induced signals in the local neigh- 
borhood or even extended portions of the Milky Way (with Milky Way 
diffusion region extending only few kpc above the Galactic plane, up 
to, possibly, 10 kpc or so [8j), spatial propagation of electrons and 
positrons is an important effect and should be appropriately modeled 
(see, e.g., 01 [II]). 



1.1.3 Convection or advection effects 

In the Galactic center region the energy loss time-scale decreases dra- 
matically, since synchrotron radiation losses are busted by the very 
large magnetic field; still, one cannot simply refer to Eq. [5j The last 
two terms on the right hand side of Eq. [2]need to be included to model, 
in the innermost part of the Galaxy, the accretion flow of gas onto the 
central black hole [H [9] . In the limit of radial free fall V is oriented 
towards the black hole and scales like — cy^ Rbh/^ (with Rbh the 
gravitational radius of the black hole, and r the radial distance) up 
to the accretion radius which is about Race ~ 0.04 pc [Jj. A particle 
propagating in such accretion flow gains momentum since it feels an 
adiabatic compression in the black hole direction. In Fig. [2] we show 
the interplay between synchrotron loss and adiabatic heating, plotting 
electron/positron trajectories in the plane radial distance versus par- 
ticle momentum, in the limit of equipartition for the magnetic field 
in the central region of the Galaxy, i.e. assuming he magnetic energy 
completely balancing the kinetic pressure jU |9] . The synchrotron loss 
dominates at high energies, while the advection gain takes over at low 
energies; electrons accumulate on the trajectory separating the two 
regimes (black curve in the figure). Since approaching the BH, the 
scaling in radius of the synchrotron loss is faster than the advection 
gain, psyn oc r~^/^ versus Padv oc r~^/^, the advection dominated re- 
gion becomes smaller and smaller and disappears for radii very close 
to the BH horizon [25]. As stated above, in the region with r > Race 
we neglect the advection and thus the trajectories are just horizontal 
lines. 



8 



10"" 10"" 10"" 10° 10" 

p [GeV] 

Figure 2: Electron/positron trajectories in the plane radius versus momen- 
tum in case of magnetic field for the Galactic center region satisfying the 
equipartition condition. Far from the turning points, synchrotron loss is 
dominant in green trajectories, while adiabatic heating takes over in violet 
trajectories. The black solid line represents the curve along which the — e~ 
accumulate since the two effects balance each other. The dotted line is the 
accretion radius Race = 0.04 pc, where advection is assumed to stop. 

A convection effect, due to the outflow of primary cosmic rays 
from the Galactic disk, where sources are located, could be relevant 
also for the propagation of particles in the whole Milky Way diffusion 
region [8j; wind velocities are expected however to be much smaller 
and their effect much less dramatic. 

1.2 Multi-frequency Emissivities and Intensi- 
ties 

The solution of Eq. [2] provides the /e~ number density in the 
stationary limit. For a radiative process i, with associated power Pj, 
the photon emissivity at frequency v is obtained by folding rig with 
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the power [27] : 



ji (x, i^) = 2 dEPi {x, E,v)ne (x , E) , (6) 

where the factor 2 takes into account electrons and positrons. The 
corresponding intensity /j, as measured by a detector, fohows from 
solution of the differential equation: 

— — — = -a{v,s)Ii{v,s,e)^ ^— , (7) 

i.e. within the increment ds along a line of sight of observation, there 
is a gain in intensity jj/(47r) ds, while a decrease ali ds could be due 
to absorption. 



1.2.1 Radio emission 

At radio frequencies, the DM-induced emission is dominated by the 
synchrotron radiation. The power for synchrotron emission takes the 
form |27] : 

Psynix, E, v) = B{x)F{v/v,) , (8) 

where mg is the electron mass, the critical synchrotron frequency is de- 
fined as Uc = 3/(4 vr) • ce/ {nieC^f B{x)E^ , and F{t) = t /~ dzK^/^iz) 
is the standard function setting the spectral behavior of synchrotron 
radiation. Synchrotron self-absorption can be a relevant effect; in- 
deed one finds that, when considering the region around the Galactic 
center, the synchrotron luminosity is so large that low frequency in- 
tensities are significantly changed by absorption [6l[71l9]; in principle, 
synchrotron self-Compton effects could be sizable as well, but they 
turn out to subdominant [9j. For all other sources, absorption can be 
safely neglected. 

It is instructive to compare the synchrotron to the gamma-ray 
luminosity under a few simplifying assumptions. Consider the limit 
of radiative losses as dominant term in the transport equation, so 
that the electron/positron equilibrium number density is given by 
Eq. [5l Suppose also that the energy loss rate is dominated by syn- 
chrotron emission, i.e. p ~ Psym and that it is sufficient to imple- 
ment the monochromatic approximation for the synchrotron power, 
i.e. assumed F{y/vc) ~ biyjv^ — 0.29) [27j. In the monochromatic 
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approximation there is a one-to-one correspondence between the en- 
ergy of the radiating electron (peak energy in the power) and the fre- 
quency of the emitted photon, that is Ep = 1/^^^(0.29 B co)~^^^ with 
Co = 3/(4 tt) • c e / (meC^)'^ , or, introducing values for numerical con- 
stants, the peak energy in GeV is Ep 2± 0.463 P'^/^S""-^/^, with 9 the 
frequency in MHz and B the magnetic field in fiG. Under this set of 
approximations, the induced synchrotron luminosity becomes: 

uLty- = J <fxp{xfEpY,{Ep) , (9) 

where we have defined Ye{E) = J^''' clE'dNe/dE' . Analogously, the 
induced 7-ray luminosity is: 

Having already stressed that electron/positron yield from WIMP an- 
nihilation is at the same level than the 7-ray yield, it follows that the 
radio and 7-ray luminosities are also expected to be roughly compa- 
rable. 



1.2.2 From the UV to the gamma-ray band 

For large magnetic field, such as for the Galactic center case, syn- 
chrotron emission may extend to the UV and (possibly) the X-ray 
band. The emission through Inverse Compton (IC) scattering of the 
ultra-relativistic electrons from WIMP annihilations on cosmic mi- 
crowave and starlight background photons give rise to a spectrum of 
photons stretching from below the extreme ultra-violet up to the soft 
gamma-ray band, peaking in the soft X-ray energy band. The Inverse 
Compton power is given by: 

/diTi 
de-^{e,x)a{e,u,E) (11) 

where e is the energy of the target photons, dn^/de is their differen- 
tial energy spectrum and a is the Klein-Nishina cross section. Addi- 
tional radiative emission is expected from the process of non-thermal 
bremsstrahlung, i.e. the emission of gamma-ray photons in the de- 
flection of the charged particles by gas electrostatic potential, as well 
as via ionization and Coulomb scattering; in general, these processes 
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Figure 3: The multi-wavelength spectrum from a 100 GeV WIMP annihilat- 
ing in bb in the dwarf spheroidal Galaxy Draco. The WIMP pair annihilation 
rate has been tuned to give a gamma-ray signal at the level of the EGRET 
limit on the gamma-ray emission. 

give subdominant contributions with respect to synchrotron and in- 
verse Compton emission. 

At gamma-ray frequencies, radiative emission matches the hard 
gamma-ray component which, as already mentioned arises from prompt 
emission in WIMP pair annihilations. We present in fig. [3] an example 
of the multi-wavelength spectrum from a 100 GeV WIMP annihilat- 
ing in bb in the dwarf spheroidal Galaxy Draco, with the details of the 
origin of radiation at various frequencies (from Ref. [13] ) . 
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2 The multi- wavelength approach and 
Galaxy Clusters 

Analyses of multi-wavelength signals from dark matter have been dis- 
cussed in the context of galaxy clusters for a few case studies, includ- 
ing the Coma cluster [11] , the Bullet cluster [12] and of the Ophiuchus 
cluster [24j. We give below a summary of the results of these studies. 

2.1 The case of Coma Cluster 

Ref. [11] performs a thorough analysis of the transport and diffu- 
sion properties of WIMP annihilation products in the Coma cluster 
of galaxies, investigating the resulting multi- frequency signals, from 
radio to gamma-ray frequencies. Other relevant astrophysical effects 
of WIMP annihilations are also discussed, including the DM-induced 
Sunyaev-Zel'dovich effect and intra-cluster gas heating. Ref. [llj finds 
that Coma radio-halo data (the spectrum and the surface brightness) 
can be nicely fitted by the WIMP-induced signal for certain particle 
physics models and under reasonable assumptions for the structure of 
the intra-cluster magnetic field. Fitting the radio data and moving to 
higher frequencies, [TT] points out that the multi-frequency spectral 
energy distributions are typically dim at EUV and X-ray frequen- 
cies (with respect to the data), but show a non- negligible gamma-ray 
emission, depending on the amplitude of the Coma magnetic field. 
A simultaneous fit to the radio, EUV and HXR data is not possible 
without violating the gamma-ray EGRET upper limit. However, the 
inverse Compton emission can reproduce both the spectrum and the 
spatial distribution of the EUV emission observed in Coma, provided 
that a quite small average magnetic field B ~ 0.15 fiG is assumed. 
Such low value of the B field is also able to make the radio data and the 
hard X-ray data of Coma consistent within a Synchrotron/IC model 
for their origins. The best-fit particle physics models yields substan- 
tial heating of the intra-cluster gas, but not sufficient energy injection 
as to explain the quenching of cooling flows in the innermost region 
of clusters. The best-fit (bb) WIMP model also yields a detectable SZ 
effect (with a peculiar spectrum very different from that of the ther- 
mal SZ effect) at the level of ~ 40 to 10 fiK in the frequency range 
~ 10 — 200 GHz, which could be observable with the next generation 
high-sensitivity bolometric arrays, space and balloon-borne microwave 
experiments, like PLANCK, OLIMPO, APEX, ALMA. 
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2.2 The case of the Bullet Cluster 



The SZuM effect is an inevitable consequence of the presence and of 
the nature of DM in large-scale structures. Ref. |12j showed that mi- 
crowave observations of the 1ES0657-556 cluster can provide crucial 
probes for the presence and for the nature of DM in cosmic structures. 
The authors calculate the expected SZ effect from Dark Matter anni- 
hilation in the main mass concentrations of the cluster, and estimate 
the sources of contamination, confusion and bias to asses its signifi- 
cance. It is found that SZ observations at w 223 GHz can resolve 
both spatially and spectrally the SZdm signal and isolate it from the 
other SZ signals, and mainly from the thermal SZ effect which is null 
at v ^ 220 — 223 GHz for the case of 1ES0657-556 . SZ observations 
of the SZjjM effect with ^ arcmin resolution and ^ /xK sensitivity 
of 1ES0657-556 are therefore crucial, and maybe unique, to find di- 
rect astrophysical probes of the existence and of the nature of Dark 
Matter, or to set strong experimental limits. 

2.3 The case of the Ophiuchus Cluster 

Ref. |24) investigates a scenario where the recently discovered non- 
thermal hard X-ray emission from the Ophiuchus cluster orig- 
inates from inverse Compton scattering of energetic electrons and 
positrons produced in WIMP dark matter pair annihilations. This 
scenario is compatible with all observational information on the clus- 
ter, and it can account for both the X-ray and the radio emission, 
provided the average magnetic field is of the order of 0.1 fiG. It is 
also shown that Fermi-LAT will conclusively test the dark matter an- 
nihilation hypothesis. Depending on the particle dark matter model, 
the LAT might even detect the monochromatic line produced by dark 
matter pair annihilation into two photons. 

3 The multi- wavelength approach and 
Dwarf Galaxies 

A class of low-background targets for multi-wavelength searches for 
WIMP dark matter annihilation is nearby dwarf galaxies. The general 
setup and an application to the specific case of the Draco dSph galaxy 
were discussed in [13] , while constraints on WIMP models from XMM 
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observations of the three satelhte dSph galaxies Fornax, Carina and 
Ursa Minor have been discussed in |21] . 

3.1 The case of Draco 

After investigating mass models for Draco in the light of available 
observational data, Ref. |13j models the dark matter density profile, 
taking advantage of numerical simulations of hierarchical structure 
formation. The gamma-ray and electron/positron yield expected for 
WIMP models is then analyzed. Unlike in larger dark matter struc- 
tures - such as galaxy clusters - spatial diffusion plays here an im- 
portant role. While Draco would appear as a point-like gamma-ray 
source, synchrotron emission from electrons and positrons produced 
by WIMP annihilations features a spatially extended structure. De- 
pending upon the cosmic ray propagation setup and the size of the 
magnetic fields, the search for a diffuse radio emission from Draco can 
be a more sensitive indirect dark matter search probe than gamma 
rays. In some cases, an extended radio emission could be detectable 
from Draco even if no gamma-ray source is identified by Fermi-LAT or 
by ACTs, making this technique the most promising search for dark 
matter signatures from the class of objects under consideration, i.e. 
nearby dwarf spheroidal galaxies. Ref. [13] also points out that avail- 
able data are consistent with the presence of a black hole at the center 
of Draco: if this is indeed the case, very significant enhancements of 
the rates for gamma rays and other emissions related to dark matter 
annihilations are expected. 

3.2 X-ray constraints from local Dwarfs 

Local dSph galaxies are an ideal environment for particle dark matter 
searches with X-rays. Ref. [21] describes how X-rays are produced 
as secondary radiation in Inverse Compton scattering off cosmic mi- 
crowave background photons of electrons and positrons resulting from 
particle dark matter annihilation. The resulting spectrum is only 
mildly dependent on the details of the particle dark matter model 
(the dark matter mass and the dominant final state into which it pair 
annihilates), and it is, generically, hard (spectral index smaller than 
^ 1.5). The normalization of the emission depends on (1) the particle 
dark matter pair annihilation rate, (2) the diffusion setup, and (3) the 
dark matter density distribution. For reasonable choices of these three 
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a priori unknown inputs of the problem, the X-ray emission is poten- 
tiahy within reach of current X-ray detectors. Interestingly enough, 
the shape of the spectral energy distribution indicates that for dSph 
galaxies X-rays have a comparable, if not better, sensitivity to indirect 
dark matter detection than gamma rays. 

Ref. |21) used XMM-Newton archival data on three Local Group 
dSph galaxies, Ursa Minor, Fornax and Carina, to search for the dif- 
fuse X-ray emission expected from dark matter annihilation. Ref. |21] 
studies the optimal energy and radial range to search for this type 
of emission, and concluded that for XMM-Newton and for the dSph 
galaxies under investigation these correspond to an energy band be- 
tween 0.5 and 8 keV and to a radius of around 6'. No significant signal 
over background is found, and this, in turn, is turned into constraints 
on particle dark matter models. The best constraints result from both 
the Fornax and the Ursa Minor observations, while data from Carina 
result in bounds that are a factor of a few weaker. Ursa Minor features 
the largest dark matter density, making it the best candidate target, 
but has the shortest usable XMM exposure. 

In determining the impact on particle dark matter searches of our 
X-ray constraints, ref. j21j points out the uncertainties resulting from 
the modeling of cosmic ray diffusion processes, and from the dark 
matter distribution. In particular, including dark matter substruc- 
tures can boost our constraints significantly. The bounds obtained 
are expressed in terms of the dominant dark matter annihilation final 
states. For those final states relevant for specific dark matter models, 
such as supersymmetry, the constraints on the mass versus pair anni- 
hilation plane are very similar. In the most conservative setup only 
models with rather large annihilation cross sections are excluded. As- 
suming a smaller diffusion coefficient, or factoring in the effect of dark 
matter substructures, our constraints fall well within the interesting 
region where the supersymmetric dark matter can be a thermal relic 
from the early universe. Also, limits can be set on particular super- 
symmetric dark matter scenarios, such as Wino or Higgsino lightest 
neutralino dark matter. 

An important result of the analysis is that even assuming a con- 
servative diffusion setup the sensitivity of X-rays and of gamma rays 
to particle dark matter annihilation in dSph galaxies are comparable. 
This fact has two-fold implications: on the one hand, if longer ob- 
servations of dSph galaxies were carried out with existing telescopes, 
it is possible that the first astronomical signature of particle dark 
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matter annihilation would come from X-rays. Secondly, should a sig- 
nature be detected in the future with gamma-ray telescopes, it would 
be extremely important to confirm the nature of the signal via X-ray 
observations. 

Future X-ray telescopes, like Constellation-X and XEUS, will also 
have greatly increased effective areas with respect to current instru- 
ments. Using the currently available projections for the effective areas 
and backgrounds of these telescopes it is estimated that the X-ray lim- 
its (0.5-8 keV band) placed by a 100 ksec observation of Ursa Minor 
with Constellation-X or XEUS would improve over the limits placed 
in this paper by factors of roughly 35 and 70, respectively. Thus even 
for a conservative diffusion model, the future generation of X-ray tele- 
scopes will place similar constraints on dark matter annihilation from 
dwarfs as Fermi-LAT, stronger constraints at particle masses below a 
few hundred GeV. In addition, a signal from Fermi could be confirmed 
with X-ray observations. 

In short. X-rays can play an important role in exploring the nature 
of particle dark matter and in pinpointing its properties. This role 
is complementary, but not subsidiary, to searches with gamma rays. 
Very exciting results at both frequencies might be just around the 
corner. 

4 The multi- wavelength approach and 
the Milky Way 

Galactic sources, including the central regions of the Galaxy as well 
as galactic DM clumps, are interesting targets for multi-wavelength 
searches for WIMP annihilations. We summarize here results on the 
innermost galactic center region from Ref. [25], dark matter "mini- 
halos" [To] and the dark matter interpretation of the so-called WMAP 
"haze" [E]. 

4.1 The Galactic Center 

Ref. [25j presents a systematic, self-consistent study of the multi- 
wavelength emission due to WIMP pair annihilations in the Galactic 
center region. High-energy electrons and positrons emitted in a region 
with large magnetic fields give rise to synchrotron emission covering 
radio frequencies up to, possibly, the X-ray band. The authors discuss 



17 



spectral and angular features, and sketch the correlations among sig- 
nals in the different energy bands. The critical assumptions for this 
region stem from uncertainties in the DM source functions, regarding 
both WIMP models and DM distributions, and from the modeling of 
propagation for electrons and positrons, and the assumptions on mag- 
netic field profiles. Radio to mm synchrotron emission is essentially 
independent from the shape of the magnetic field in the innermost 
region of the Galaxy, while at shorter wavelengths, i.e. in the infrared 
and, especially, the X-ray band, a different choice for the magnetic 
field may change predictions dramatically. Radio signals have in gen- 
eral very large angular sizes, larger than the typical size for the source 
function and hence of the 7-ray signals. The size of the region of 
synchrotron X-ray emissivity shrinks dramatically going to larger fre- 
quencies, smaller WIMP masses, or softer annihilation channels. 

The luminosity of the WIMP source at the different frequencies, 
and especially comparing the radio to the 7-ray band, is essentially 
at a comparable level, with luminosity ratios depending rather weakly 
on WIMP mass and annihilation channels. This is interesting, since 
the GC astrophysical source Sgr A*, an unusual source, certainly very 
different from typical galactic or extragalactic compact sources associ- 
ated to black holes, has a very low luminosity over the whole spectrum, 
at a level at which it is plausible that a WIMP-induced component 
may be relevant. Indeed, while none of the fluxes detected in GC di- 
rection has spectral or angular features typical of a DM source, still 
all data-sets contribute to place significant constraints on the WIMP 
parameter space. Although the 7-ray band is the regime in which 
it is most straightforward to make the connection between a given 
dark matter model and the induced signal (hence it is also the regime 
on which most of previous analyses have concentrated on), it does not 
seem to be the energy range with the best signal to background ratios. 
In the case of large magnetic fields close to the GC, X-ray data can 
give much tighter constraints. Radio and NIR measurements, which 
are less model dependent, tend to be more constraining as well. 

4.2 Galactic Mini-halos 

Ref. [10] calculates the spectrum and spatial extent of diffuse emis- 
sion from the charged particle products of dark matter annihilations in 
galactic satellites that are currently within the diffusion zone, namely 
within a few kpc of the stellar disk. For satellites moving with typical 
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galactic halo velocities of 300 km si , the crossing time of the diffu- 
sion zone is of the same order as the diffusion time, thus an inherently 
time-dependent treatment is required. For annihilation sources, e.g. 
galactic satellites at typical distances of 10 kpc, the diffuse emission 
in both inverse Compton and synchrotron extends over roughly 300 
square degrees. At least in terms of the number of photons, the dif- 
fuse inverse Compton emission might be detectable by Fermi-LAT, 
assuming bright enough annihilation sources. The spatial extent of 
the emission makes however its detection problematic. 

4.3 The WMAP "Haze" 

Ref. |15j shows that the observed features of the WMAP haze match 
the expected signal produced through the synchrotron emission of dark 
matter annihilation products for a model with a cusped halo profile 
scaling as p(r) oc r^^'^ in the inner kpc, and an annihilation cross 
section of ~ 3 x 10"^^ cm^/s. A wide range of annihilation modes are 
consistent with the synchrotron spectrum, and no boost factors are 
required. The properties required of a WIMP to generate the haze 
are precisely those anticipated for the most theoretically attractive 
particle dark matter candidates. If the haze is generated through 
dark matter annihilations, this will have very interesting implications 
for Fermi. If the p{r) cx r~^'^ slope of the halo profile continues to 
the inner Galaxy, the gamma ray flux from the Galactic Center will 
be observable by Fermi, so long as the WIMP is lighter than several 
hundred GeV, in spite of the presence of the observed HESS source in 
the region. 

5 Radio Observations 

We summarize here the findings of two studies, [20] and [26], that 
addressed the detection of a radio emission from WIMP annihilation. 

Ref. |20] calculates the secondary WIMP-induced radio emission 
from the galactic halo as well as from its expected substructures, com- 
paring it with the measured diffuse radio background. They employ a 
multi-frequency approach using data in the relevant frequency range 
100 MHz-100 GHz, as weh as the WMAP Haze data at 23 GHz. The 
derived constraints are of the order {av) = 10~^^ cm^s~^ for a DM 
mass = 100 GeV sensibly depending however on the astrophysical 



19 



uncertainties, in particular on the assumption on the galactic mag- 
netic field model. The signal from single bright clumps offers only 
weak signals because of diffusion effects which spread the electrons 
over large areas diluting the radio signal. The diffuse signal from the 
halo and the unresolved clumps is instead relevant and can be com- 
pared to the radio astrophysical background to derive constraints on 
the DM mass and annihilation cross section. Conservative constraints 
are at the level of {av) ~ 10~^^ cm^s~^ for a DM mass = 100 GeV 
from the WMAP Haze at 23 GHz. However, depending on the astro- 
physical uncertainties, in particular on the assumption on the galactic 
magnetic field model, constraints as strong as {av) ~ 10"'^^ cm'^s"^ 
can be achieved. 

Ref. [26] calculates the intensity and angular power spectrum of the 
cosmological background of synchrotron emission from WIMP anni- 
hilations into electron positron pairs, and compares this background 
with intensity and anisotropy of astrophysical and cosmological ra- 
dio backgrounds, such as from normal galaxies, radio-galaxies, galaxy 
cluster accretion shocks, the cosmic microwave background and with 
Galactic foregrounds. Under modest assumptions for the dark matter 
clustering, Ref. [26] finds that around 2 GHz average intensity and 
fluctuations of the radio background at sub-degree scales allows to 
probe dark matter masses ^ 100 GeV and annihilation cross sections 
not far from the natural values {av) ~ 3 x 10~^^ cm^ s~^. The angular 
power spectrum of the signal from dark matter annihilation tends to 
be flatter than that from astrophysical radio backgrounds. Further- 
more, radio source counts have comparable constraining power. Such 
signatures are interesting especially for future radio detectors such as 
SKA. 

6 Conclusions and Overview 

If dark matter particles pair annihilate into stable Standard Model 
particles, a population of energetic, non-thermal electrons and positrons 
is necessarily produced. The secondary radiation resulting from the 
energy losses of this non-standard population (including synchrotron 
and inverse Compton up-scattering of background radiation photons) 
in turn produces a peculiar multi-wavelength spectrum, extending 
from radio to gamma-ray frequencies. In some cases, constraints 
from radio or X-ray observations of galactic or extragalactic targets 
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place constraints on dark matter models which are more stringent 
than those resulting from gamma-ray data. Reversing this statement, 
multi-wavelength observations will complement gamma-ray observa- 
tions as probes of particle dark matter, since the expected luminosi- 
ties at different frequencies are generically comparable. The indirect 
search for dark matter with astronomical observations at various fre- 
quencies is therefore a crucial and indispensable element in the quest 
for the fundamental nature of dark matter. 
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